Isotropy theorem for cosmological vector fields 



(N 

o 

(N 

\o 

O 
U 

Oh: 

I 

o 

(N 

> 

(N 

\q 

rn 

O 



X 



J. A. R. Ccmbranos, C. Hallabrin, A. L. Maroto and S. J. Nunez Jarcno 
Departamento de Fisica Teorica I, Universidad Complutense de Madrid, E-28040 Madrid, Spain 

(Dated: March 1, 2013) 

We consider homogeneous abehan vector fields in an expanding universe. We find a mechanical 
analogy in which the system behaves as a particle moving in three dimensions under the action of a 
central potential. In the case of bounded and rapid evolution compared to the rate of expansion, we 
show by making use of the virial theorem that for arbitrary potential and polarization pattern, the 
average energy-momentum tensor is always diagonal and isotropic despite the intrinsic anisotropic 
evolution of the vector field. For simple power-law potentials of the form V = \{A^ Afi)" , the 
average equation of state is found to be ui = (n — l)/(n + 1). This implies that vector coherent 
oscillations could act as natural dark matter or dark energy candidates. Finally, we show that under 
very general conditions, the average energy-momentum tensor of a rapidly evolving bounded vector 
field in any background geometry is always isotropic and has the perfect fiuid form for any locally 
inertial observer. 



Our knowledge about the history of our universe has 
improved over the last years with the advent of a large 
amount of new observations. There are robust astro- 
physical data that support the existence of an early in- 
flationary era; an additional matter component supple- 
menting the baryonic one, known as dark matter; and 
a present era of accelerated expansion driven by dark 
energy. However, the fundamental nature of these com- 
ponents remains unknown. Different phenomenological 
and model building strategies have been proposed to un- 
derstand them in which rapidly evolving fields could play 
an important role. 

In particular, different scalar dynamics have been con- 
sidered as the possible answers for these open questions: 
standard inflation models arc based on the slow-roll evo- 
lution of a scalar field known as inflaton. The dynam- 
ics of this field is able to finish successfully the accel- 
erated regime and its later rapid oscillations provide a 
mechanism for reheating, transferring perturbatively or 
non-perturbatively the energy density of the oscillations 
to the matter fields. Also the possibility of generating 
inflationary expansion during the period of inflaton os- 
cillations has been studied in [l|. Another example of 
rapidly evolving scalar fields can be found in non-thermal 
dark matter candidates like the axion [2|1 or other mas- 
sive scalar [3| or pseudoscalar particles [j|, in which the 
energy density of rapid scalar coherent oscillations scales 
precisely as non-relativistic matter Q . Dark energy mod- 
els based on the dynamics of scalar fields are commonly 
known as quintessence. Oscillating evolutions have also 
been considered within this context [6|. 

On the other hand, it is interesting to remark that all 
these possibilities could in principle be offered by any 
bosonic degree of freedom and not only by scalars. In 
fact, a large number of fundamental vector fields arc 
present in the standard model of particles and interac- 
tions, and in the most part of its extensions. Therefore 
it is natural to consider vector models which could shed 
light on the above mentioned open problems in cosmol- 
ogy. However, there is an important distinctive feature 
of vector fields as compared with scalar fields. Even in a 



homogeneous configuration, vector dynamics is generally 
anisotropic. 

Most part of observational data seem to be consistent 
with an early homogenous and isotropic universe. Mod- 
els supporting a large amount of anisotropy suffer severe 
constraints. However, despite this fact, there are some 
examples in which vector fields have been shown to pro- 
vide interesting models in different cosmological scenar- 
ios. For instance, models of inflation based on vector 
fields have been extensively studied recently [3], and a 
first proposal for vector inflation which could avoid the 
generation of an excess of anisotropy can be found in [8| . 
Models of dark energy based on massive vector fields have 
been also considered in Q. Vector dark energy without 
potential terms has been proposed in jlOj . Vector mod- 
els for dark matter based on hidden sector gauge bosons 
have been discussed in [11[. A possible role in the gen- 
eration of metric perturbations in the so called curvaton 
scenario, has also been considered for vectors in |12| . In 
that work, it was shown that an oscillating massive ho- 
mogeneous vector field behaves as non-relativistic matter 
with an equation of state w = 0, in a completely anal- 
ogous way to the scalar case. Despite the anisotropy 
of the oscillations, the average energy-momentum tensor 
(EMT) turned out to be isotropic in that case. This fact 
has been used in order to propose oscillating massive vec- 
tor fields as non-thermal dark matter candidates in |13| . 

On the other hand, the potential of anisotropic cos- 
mologies have received important attention in the last 
years, mainly motivated by the possible existence of 
anomalies in the isotropy of the Cosmic Microwave Back- 
ground (CMB) and matter distributions which could be 
pointing to the existence of a preferred spatial direction 
in the Universe |14| . 

In this work, we show that despite its intrinsically 
anisotropic evolution, the average EMT associated to 
rapidly evolving vector fields is isotropic under very gen- 
eral and natural conditions. The proof has a mechanical 
analogy with the virial theorem, and applies to a general 
Abclian gauge vector with any general self-interaction 
given by a potential of the form V{A^A^^). 



In order to simplify the argument, let us consider first 
a flat Friedmann-Lemaitre-Robertson- Walker (FLRW) 
metric given by: 



ds^ = dr -a\t)dx 



(1) 



The Lagrangian density for a vector field with a potential 
which is an arbitrary scalar function of A^ = A^A^^ is 
given by: 



C = -^F^,F^''-V{A^) 



(2) 



where, for an Abelian field, the field strength tensor is 

F^. = a^A, - d.Af, . (3) 

The corresponding field equations read: 

i^^% + 2V'{A^)A'' = , (4) 

where V'{x) = dV/dx. We also need to calculate the 
EMT: 

+ V{A^)g^^ ~ 2V\A^)A^'A^ . (5) 

Assuming homogeneity for the vector field, i.e. A^ = 
(ylo(t), j4i(i)), we obtain for the temporal component: 

V'{A^)A^ = , (6) 

whereas for the spatial components, we get: 

A, + HA,~2V'{A^)A,^Q. (7) 

The EMT components read: 






2 a2 



(8) 



^^^-^^-2 a^ 



S^j 



A,A 



k^k 



ViA^) - 2V'iA^)^^, fc = 1, 2, 3 ; (9) 



T^o = ; 



rjit 



J\i-/\.j 



+ 2V'{A') 



11 A^\ '^^ 



i + 3 



(10) 
(11) 



Notice that in the definition of the pressures pfe, no sum- 
mation in k is assumed. 

Now we can express the conservation law T^^.^ = as: 



p + i7K^Pfc + 3p =0 



(12) 



The temporal equation ([61), implies"'^ 



(13) 



Therefore the homogeneity condition reduces the prob- 
lem to the evolution of a 3- vector A[i). We can define 
the typical time evolution scale of the spatial compo- 
nent Ai from: Wi ~ |Ai/Ai|. This value coincides with 
the frequency for a oscillatory movement. In the case of 
rapid evolution of Ai in relation to the universe expan- 
sion uji':^ H ,we can neglect time derivatives of the scale 
factor so that defining r^ = Ai/a in a time range of order 



we have: 



A, 



Hr, 



a 



(14) 



Thus, we can ignore the friction term and we see that Eq. 
([7|) reduces to the evolution equation of a point particle 
with position vector f = A/a in the presence of a cen- 
tral potential V{—r'^). In particular, this implies that we 
can make use of the standard classical mechanics results. 
Thus, since the potential is central, we will have conserva- 
tion of the corresponding angular momentum i = r x r, 
which, in turn, implies that the vector A{t) should evolve 
in a fixed plane, orthogonal to L. 

On the other hand, in the case of rapid evolution we 
can also write: 



Thus, the total energy density can be written as: 



(15) 



P=l^S^^+ViA' 
2 a 



¥ 



^ + V{-r^), (16) 



which can be considered as constant within the short 
time-scale of variation of the vector field. 

Thus, we are left with the reduced radial problem, i.e. 
the motion of a particle in the radial dimension in the 
presence of the effective potential: 



Veffir) 



2r2 



+ V{-r^) 



(17) 



where the value of the constant L is set by the initial 
conditions. The other constant of motion associated to 
this problem is the energy density po = Vm- 

As commented before, if the temporal evolution of Ai 
is sufficiently rapid and we concentrate only in a time 
interval of order uJi~^, we can ignore the effect of the 
universe expansion, i.e. we can consider a as constant. 



We see that the off-diagonal part of the EMT does not 
contribute in (|12p because of the homogeneity of the vec- 
tor field. 



^ There is another possible solution: V'(A'^) = 0. In this case, 
the isotropy theorem cannot be applied because the evolution 
of Ai is not rapid. It can be showed that Ai = Cj/a, and the 
anisotropic components of the EMT decay as a"* [15|| . 



In this case, we can write the equation of motion for the 
vector field as Ai = 2V'{A'^)Ai. Let us define: 



G 



AiAj 



«,J = 1,2,3; 



and take its time derivative: 
G. 



Ai Aj Ai Aj 



''" a^ 



Using the equations of motion, we get: 

. AjAi AiA 



G,, - 2V'{A') 



I { a'^\ ^ 1 



i^j 



(18) 



(19) 



(20) 



Integrating this expression in a given time interval [0, T] 

AiAj 



Gij{T)-G.M 



1 1 /i2\ ^«^J 



2V\A^) 



/x-i-ixj 



(21) 



with i = 1,2,3. If the motion is periodic and T corre- 
sponds to the oscillation period, the left hand side (l.h.s.) 
vanishes. If the motion is not periodic, but Ai and Ai 
are bounded, by taking T sufficiently large, but satisfy- 
ing H~^ ^ T ~:^ (jj~^ for any i, we can also neglect the 
l.h.s. of the equation. In those cases, we have: 



^) = -(2F'(A2)^), *,j = 1,2,3. (22) 



By using this result in (|11|) . we straightforwardly get: 

(t;.)=0, i^j- (23) 

i.e. the average EMT is diagonal. Also, using ^F^ in ([5]), 
we obtain 

1 AiAi 



{p,)^-{Tl) = lt^5^,\^(^v{A^)) 



2 a2 



k = 1,2,3 



(24) 



i.e. all the average pressures are the same, so that the 
isotropy of the mean value of the EMT is proved: 



(25) 



{T)) = {p) 5\ 



It is interesting to note that the virial relations 
also allow us to get the average equation of state for a 
power-law potential: V = A(A^A^)". In this case, we 
have: 



1 AjAj_ y 

2 a2 



\V{A^)) 



(26) 



what implies: 

(p,) = (p) = (n - 1) (y(A2)) , fc = l,2,3. (27) 
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FIG. 1: The effective equation of state a; = (p) / (p), can be 
computed numerically by integrating expressions ([8]) and (J9l) 
with respect to tlie r variable with the use of Eq. (|15p and 
P6p . Thus, in the case V — —aA^, uj — and the coher- 
ent vector oscillations behave as cold dark matter (CDM); or 
for V = 6yl*, oj — 1/3 and they behave as radiation. These 
results agree with Eq. (|29|) . For a general potential, the ef- 
fective equation of state depends on the initial conditions: 
Vm and L. In this figure, one can see the L dependence of 
u) ior V = -aA^ + M* and Vm = lO". Depending on the 
particular values of a and b, the equation of state interpo- 
lates between the radiation and the CDM behavior; from 
top to bottom respectively, dotted (black) line corresponds 
to (a, 6) = (10~^,2.4 • 10~^), blue (dashed) line corresponds 
to (a,fe) = (0.1,2 ■ 10""), dashed-dotted (red) for (a, &) = 
(0.4,10"''), long dashed (orange) for (a,fe) = (0.7,6 ■ 10"^) 
and continuous (green) for (a, b) = (10"^, lO"'^^). 



On the other hand, for the average energy-density, we get 
from m-. 



{p)^{n+l){V{A^)) 



and finally: 



M 



1 ' 



(28) 



(29) 



which agrees with the scalar case Q and shows that the 
average equation of state does not depend on the polar- 
ization of the vector oscillations (notice that the scalar 
case corresponds to L = 0, i.e the evolution of a particle 
in one spatial dimension). In the case of a general poten- 
tial, not necessarily of the power-law form, the equation 
of state cannot be obtained analytically and, in general, 
it will depend on the initial conditions for the vector oscil- 
lations, as Fig. 1 shows. Notice that the average energy 
density and pressure satisfy the conservation equation 
(fT2|) up to corrections of order 0{HT). 

We could think that the isotropy of the average EMT 
could have been inherited from the isotropy of the FLRW 



metric. But repeating the process for a Bianehi I metric: 

ds^ = dt^ - al{t)dx^ - al{t)dy^ - al{t)dz^ , (30) 



if we replace 



^i -^i 



-, 1 = 1,2,3 



(31) 



and assume w^ much greater than Hi ~ bujai^ it is 
straightforward to obtain the same results for (T^^). 

The previous results can be directly extended to 
general space-time geometries (not necessarily homoge- 
neous). With that purpose, let us consider a locally in- 
ertial observer at Xq = and write the metric tensor in 
Riemann normal coordinates around Xn 11611: 



Let assume that the following conditions hold: 



(32) 



1. The Lagrangian of the vector field is restricted to 
be given by the form of Eq. ([2]) . 

2. The vector field evolves rapidly: 
\Rl^.^<■^l and \d,A,\ <C |i.|, 1,3 = 1,2,3 (33) 

for any component of the Riemann tensor. 

3. Ai and Ai remain bounded in the evolution. 

The second condition implies that if we are only inter- 
ested in time scales of order ujj , then we are in a normal 



neighborhood and we can neglect the second term in ((32)) 
and also work with a homogeneous vector field. In such 
a case, it is possible to rewrite all the above equations 
with a = 1. Thus, by using an interval [0, T] that verifies 
the condition: 



\Rl,.\ 



<r 



-^«a;? 



(34) 



for any of the components of the Riemann tensor and 
any of the spatial components of the vector: z = 1, 2 and 
3, it is possible to obtain (PT|) and prove that the mean 
value of the EMT is isotropic. It is interesting to note 
that if the metric is non-homogeneous, the EMT can be 
non-homogeneous, but its average value will be isotropic 
as seen from a locally inertial frame. 

In this work, we have considered the evolution of gen- 
eral Abelian vector fields in an expanding universe. We 
have shown by means of a mechanical analogy with a par- 
ticle moving in a central potential that a generalized ver- 
sion of the virial theorem ensures that, for rapid dynam- 
ics, the average EMT is always diagonal and isotropic, 
despite the fact that the evolution always takes place in 
a fixed plane. The result can be extended to arbitrary 
geometries within a normal neighborhood of any locally 
inertial observer. 
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